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We use the newly published 28 observational Hubble parameter data (H(z)) and current largest 
SNe la samples (Union2.1) to test whether the universe is transparent. Three cosmological- 
model-independent methods (nearby SNe la method, interpolation method and smoothing method) 
are proposed through comparing opacity-free distance modulus from Hubble parameter data and 
opacity-dependent distance modulus from SNe la . Two parameterizations, r(z) = 2ez and 
. t(z) = (1 + z) 2e — 1 are adopted for the optical depth associated to the cosmic absorption. We 

, find that the results are not sensitive to the methods and parameterizations. Our results support a 

ON . transparent universe. 
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O 1. INTRODUCTION 

O 

The unexpected dimming of Type la supernova (SNe la) is thought to be the evidence of acceleration of the universe 
In the frame of General Relativity (GR), the most famous explanation is the existence of dark energy with a 
negative pressure However, there are some issues on this plausible mechanism for observed SNe la dimming. The 
. photon number conservation may be deviated. For example, it is due to the dust in our galaxy and oscillation of 

(N ■ 

I i photons propagating in extragalactic magnetic fields into very light axions. These absorption, scattering or axion- 
photon mixing may lead to dimming [jj. Other mechanisms are widely proposed including modified gravity Q], 
Q] , evolutionary effects in SNe la events Q] , violation of cosmological principle Q 
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^vq dissipative processes |5j , evolutionary effects in SNe la events [6| , violation of cosmological principle [7( and so on. On 
the other hand, the deviation of photon number conservation is related to the correction of Tolman test which is 
equivalent to measurements of the well-known distance-duality relation (DDR) 



where Dr, is luminosity distance, Da is angular diameter distance and z is redshift. The DDR is in fact a particu- 
lar form of presenting the general theorem proved by Etherington known as the "reciprocity law" or "Etherington 
reciprocity theorem" . DDR holds for general metric theories of gravity in any cosmic background and it is valid for 
any cosmological models based on the Riemannian geometry. It is independent of gravity equation and the universe 
components. However, DDR may be not valid in the case that photons do not travel along null geodesies or the 
cosmic opacity exists. Many efforts have been done to test DDR though astronomical observations [10j. Usually, they 
assume the form 75~(1 + z ) 2 — 7 l{ z ) i where rj{z) = 1 + t]qz or rj(z) = 1 + ?/oy^I- Compared to conservation of 



* liaokai@mail.bnu.edu.cn 



2 



photon number, the assumptions that the mathematical tool used to describe the space-time of universe is Riemannian 
geometry and photon travels along null geodesic are more fundamental and unassailable, thus the deviation of DDR 
most possibly indicates cosmic absorption. In this case, the flux received by the observer will be reduced by a factor 
e~ T ( z \ and observed luminosity distance can be obtained by 11 [ 

D L ,obs = D L . true e T ^' 2 , (2) 

where r(z) is the optical depth related to the cosmic absorption. The relation between t(z) and r/(z) is e r ( z )/ 2 = r](z) 



12 1 . Following this assumption, Avgoustidis et al. [13] studied the difference between SNe la observations and Hubble 
parameter data. H(z) data are mainly obtained through the measurements of differential ages of red-envelope galaxies 
known as " differential age method" . The aging of stars can be regarded as an indicator of the aging of the universe. 
The spectra of stars can be converted to the information of their ages, as the evolutions of stars are well known. 
Since the stars cannot be observed one by one at cosmological scales, people usually take the spectra of galaxies 
which contain relatively uniform star population. Moreover, H(z) data can be obtained from the BAO scale as a 
standard ruler in the radial direction known as " Peak Method" . These methods are apparently independent of galaxy 
luminosity so that it will not be affected by cosmic opacity. However, SNe la observations are affected by many 
sources of opacity, such as the hosting galaxy, intervening galaxies, intergalactic medium, the Milky Way and exotic 
physics which affect photon conservation. Under the assumption Dl — Aa(1 + z ) 2+e , they investigated the cosmic 
opacity by confronting the standard luminosity distance in spatially flat ACDM model with the observed one from 
SNe la observations. Combining with the H(z) data, which is not affected by transparency but yields constraints 
on f2 m , and marginalizing over all other parameters except e, they got e = — 0.04 j^g® (2cr). Noticing that their 



method depends on cosmological models, Holanda et al. 14{ further proposed a model- independent estimate of Dl 
which are obtained from a numerical integration of H(z) data. They also explored the influence of different SNe la 
light-curve fitters (SALT2 and MLCS2K2) and found a significant conflict. Based on Holanda et al., we proposed 
three model-independent methods that are different from theirs to explore the cosmic opacity. They firstly got the 
luminosity distances from H(z) data (12 data) at corresponding redshifts (at H(z) data) and gave a polynomial fit 
based on these 12 luminosity distances data with their errors, then calculated the values at the redshifts corresponding 
to SNe la through this polynomial fit curve. On the contrary, we use SNe la data to get the luminosity distances 
at the redshifts corresponding to H(z) data through interpolation method, smoothing method and nearby SNe la 
method. Our data sets contains 28 available H(z) data and the largest SNe la samples Union2.1 15 1. 

The Letter is organized as follows: In Section 2, we introduce the method of obtaining luminosity distance from 
H(z) data. In Section 3, we give the three methods that can convert SNe la luminosity distances to the luminosity 
distances at the redshifts corresponding to H(z) data. In Section 4, the results are performed. Finally, we make a 
conclusion in Section 5. 



2. LUMINOSITY DISTANCE FROM OBSERVATIONAL HUBBLE PARAMETER DATA 



In this section, we introduce the method proposed by Holanda et al. 
can be written in this form 

*{*) = - 1 ^ 



1 + z dt 



14j . The expression of the Hubble parameter 

(3) 



17J used the 



which depends on the differential age as a function of redshift. Based on Jimenez et al. [if]], Simon et al. 
age of evolving galaxies and got nine H(z) data. Stern et al. [3] revised these data at 11 redshifts from the differential 
ages of red-envelope galaxies. Gaztanaga et al. 1131 took the BAO scale as a standard ruler in the radial direction, 
obtained two data. Recently, Moresco et al. [20j obtained 8 data from the differential spectroscopic evolution of 



early-type galaxies as a function of redshift. Blake et al. 2l| obtained 3 data through combining measurements of the 
baryon acoustic peak and Alcock-Paczynski distortion from galaxy clustering in the WiggleZ Dark Energy Survey. 
Zhang et al. [22] obtained another 4 data. Totally, we have 28 available data summarized in Table 1. Following 



Holanda et al. 



14j . we transform these 28 H(z) data into luminosity distance. Using a usual simple trapezoidal rule, 



the comoving distance can be calculated by 
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Since we use much more data than Holanda et al. (they used only 12 data), this trapezoidal rule will work much better. 
In Fig. 1, we show the relative error with respect to the data number at a characteristic redshift z = 1. We assume 
a standard ACDM model with tt m — 0.3, J7a = 0.7, and divide z — 1 into different numbers of intervals averagely, 
then calculate the relative errors according to Eq. (4). We find that the relative errors decrease remarkably when the 
numbers of intervals increase from 12 to 28. With the standard error propagation formula, the error associated to the 



bin is given by 



1/2 
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(5) 



where a Hi is the error of H(z) data. The error corresponding to certain redshift is the sum of Sj. The Hubble constant 
Ho = 73.8 ± 2.4 km/s/Mpc 23j is used in our study. The 28 luminosity distance data from H(z) data are shown in 
Fig. 2, as well as the Dl from Union2.1 SNe la samples. 
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FIG. 1: Relative error with respect to the data number at z = 1. tt m = 0.3, Qa = 0.7. 
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FIG. 2: Dl(z) obtained from measurements of the Hubble parameter data and the Union2.1 samples, respectively. 

3. DEALING WITH SNE IA SAMPLES 

In this section, we introduce three methods through which we can obtain the luminosity distance of one certain 
SNe la point at the same redshift of the corresponding H(z) data. 

3.1. nearby SNe la method 

Since the SNe la Samples is much larger than H(z) data, the nearby SNe la can be substituted for the one at the 
redshift corresponding to H(z) data. Points zsNe ia — zh are centered around the line Az = 0, as shown in Fig. 
3 which plots the subtractions of redshifts between H(z) data and the associated SNe Ia. Similar with the DDR 
test we have to choose a criterion based on the data. Our selection criterion is Az = \zh — ZsNe ia\ < 0.003. 
This selection criteria can be satisfied for most of the H(z) data except for the points at z = 0.9 and z = 1.037 
(z = 1.43,1.53,1.75 are obviously ruled out) and small enough to reduce the systematic errors and guarantee the 
accuracy. 




3.2. interpolation method 

In order to avoid any bias brought by redshift incoincidence between H(z) data and SNe Ia, as well as to ensure 
the integrity of the H(z) data, we can use the nearby SNe Ia points to obtain the luminosity distance of SNe Ia 
point at the same redshift of the corresponding H(z) data. This situation is similar with the cosmology-independent 
calibration of GRB relations directly from SNe Ia 
and the error can be calculated by 



24j. When the linear interpolation is used, the distance modulus 



(i = fjn + (fx i+1 - m)[{z - Zi)/(z l+1 - Zi)]. 



(6) 
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FIG. 3: Subtractions of redshifts between -ff(z) data and the associated SNe la. 



<r M = ([(z 4+1 - z)/(z i+1 - Zi )f^ + [(z - Zi )/(z i+1 - ^)] 2 < l+1 ) 1/2 , (7) 

where subscripts i and i + 1 stand for the nearby data. We use the same method as used in [25] to obtain the best 
estimate for each SNe la which is weighted average of all available distance moduli at the same redshift. The derived 
distance modulus for each SNe la is 



m = (Ew/0/(EO> (8) 

i i 

with its uncertainty = 0W 2 ) -1 ' 2 - 



3.3. smoothing method 



We introduce a non-parametric method of smoothing supernova data over redshift using a Gaussian kernel in 
order to reconstruct luminosity distance [261]. This procedure was initially used in the analysis of the cosmic large 



scale structure 



27] . Through this mo del- independent method, we can extract information of various cosmological 



parameters, such as Hubble parameter, the dark energy equation of state, the matter density. Wu and Yu [28 1 
generalized this method to eliminate the impact of Hq and obtained the evolutionary curve of luminosity distance 
using SNe la Constitution set and Union2 set. In this Letter, we follow this generalized method to get the luminosity 
distance curve using Union2.1 set. Firstly, we obtain the variable ln/(z) = lnD L (z) — \nh through iterative method 



\nf(z) s n = In f{ Z )' n _ x +N{z) V [\nf° bs ( Zl ) - In /(z)^] x 



exp 
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where the reduced Hubble constant h = i^o/100 and the normalization factor 

N(z)- 1 =$>xp 
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(10) 



The value of parameter A was discussed by Shafieloo et al. [26j . The results are not sensitive to A. Following 



Shafieloo et al. 



261 ] and Wu and Yu 28j, A = 0.6 is used in this Letter. Eq. (9) can give the smoothed luminosity 



distance at any redshift z after n iterations. When n = 1 



hi f(z)l= In f(zy Q + N(z)J2[^f obs ( Zl )- In f(z) s ] x exp 



= In D L (z) s + N(z) J2 [ln/° bs (^) - In D L (z) s ] x exp 



In 2 ( TT^ 
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(11) 



whereDi(2;)g depends on the suggested background cosmological model. Following Wu an Yu [28J, we adopt uCDM 
model with u> = —0.9 and fl m o = 0-28 as the background model. The relation between \nf obs (zi) and h\D° ] bs {zi) 
which is obtained from observed SNe la is 



ln/ obs (^) 



In 10 



[^ bs { Zi ) - 42.38] = \nD° L bs (z t ) - In h, 



(12) 



A* (zi) is the observed distance modulus from SNe la. In order to know whether we get a best-fit value after some 
iterations, we calculate, after each iteration, % 2 



2 =Sr (j*(z i ) n -tl° bs (z i )) 2 



a 



(13) 



The best- fit result is corresponding to the minimum value of X s n- The 1 a corresponds to A% 2 = 1. The results For 
Union2.1 samples are shown in Fig. 4 and Fig. 5, the minimum value n = 32. 



CONSTRAINTS ON COSMIC OPACITY 



The observed distance modulus can be expressed as 11] 



Vobs{z) = fitrue(z) + 2.5[log 10 e]r(z) . (14) 

To examine the sensitivity of test results on the parametric form, we adopt two parameterizations, r(z) — 2ez and 
t(z) = (1 + z) 2t — 1 which are not strongly wavelength dependent on the optical band [29j . e here describes the cosmic 
opacity. The former one is linear and it can be derived from the DDR parameterization = Da(1 + z) 2+t for small 
e and redshift. To constrain the value of e, we use the usual maximum likelihood method of \ 2 fitting 

2 _ \ " iHobsji) - utrueji) - 2.5[log 10 e]r(^)) 2 
x ^77T+^ \ ' 

I /j,{obs) fi(true) 
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FIG. 4: Computed Xs f° r the reconstructed results at each iteration. 
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FIG. 5: Smoothed luminosity distance for Union2.1 samples. 

the subscript i stands for the data at the redshifts corresponding to H (z) data. Our results are shown in Fig. 6, Fig. 
7 and Fig. 8, respectively. For smoothing method, we consider two cases: i m ax = 25 and i m ax — 28 (containing the 
data at z = 1.43, 1.53, 1.75). From the likelihood of e using different methods, we can see current SNe la samples and 
H(z) data support a transparent universe. These results are slightly different from Holanda et al. [3] while their 
results seem a little prone to a non-transparent universe especially with MLCS2K2 compilation. 
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FIG. 7: Likelihood function for e using interpolation method. 



5. CONCLUSION 



Until now, modern cosmology has discovered many interesting phenomena behind which there exists underlying 
physical mechanisms. In principle, since the current astronomical observations are not precise enough to distinguish 
between different cosmological models, for example, various dark energy models are consistent with observations, we 
can explore all the possibilities among which the true one exists. Though the matter component in the universe is so 
diluted, photons will get though a huge space to observers, photon conservation can be violated by simple astrophysical 
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FIG. 8: Likelihood function for e using smoothing method. 



effects or by exotic physics. Amongst the former, the attenuation is due to interstellar dust, gas, plasmas and so on. 
More exotic sources of photon conservation violation involve a coupling of photons to particles beyond the standard 
model of particle physics. Therefore, the concept of cosmic opacity should be considered naturally. In this Letter, we 
use the current observational Hubble parameter data which is opacity-free and SNe la observations which depends 
on cosmic opacity to test whether the universe is transparent. The results from three model- independent methods 
converge to a point that the effects of cosmic opacity is vanished. For future study on this problem, we think 
the wavelength is a considerable factor and more independent methods of testing cosmic opacity will confirm the 
conclusion. 
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